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Abstract 

We develop a direct method of solution for finding the bright iV-soliton solution of the 
Fokas-Lenells derivative nonlinear Schrodinger equation. The construction of the solution 
is performed by means of a purely algebraic procedure using an elementary theory of 
determinants and does not rely on the inverse scattering transform method. We present 
two different expressions of the solution both of which are expressed as a ratio of determi- 
nants. We then investigate the properties of the solutions and find several new features. 
Specifically, we derive the formula for the phase shift caused by the collisions of bright 
solitons. 
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1. Introduction 

In this study, we consider the following Fokas-Lenells (FL) derivative nonlinear Schrodinger 
(DNLS) equation: 

u xt = u - 2i\u\ 2 u x , (1.1) 

where u = u(x, t) is a complex- valued function of x and t, and subscripts x and t appended 
to u denote partial differentiations. The original version of the FL DNLS equation has 
been derived as an integrable generalization of the NLS equation using bi-Hamiltonian 
methods [1]. The equation (1.1) is a gauge transformed form of it. We employ the above 
simplified equation for convenience. Recently, the equation was shown to model the 
propagation of nonlinear light pulses in monomode optical fibers when certain nonlinear 
effects are taken into account [2]. It admits a Lax pair which enables one to solve the initial 
value problem by means of the inverse scattering transform (1ST) method [3]. Actually, 
a few conservation laws and a bright 1-soliton solution were obtained explicitly in [3]. 
Subsequently, the general bright iV-soliton solution was constracted with the use of the 
dressing method [4], where N is an arbitrary positive integer. Note, however, that the 
detailed analysis of the multisoliton solutions has not been undertaken. For the purpose 
of obtaining soliton solutions in a simple manner, one can employ a direct approach which 
is now known as Hirota's direct method or bilinear transformation method [5, 6]. The 
method has been applied successfully to a large number of soliton equations including the 
Korteweg-de Vries, NLS and sine-Gordon equations. 

The purpose of the present paper is to construct the bright iV-soliton of the FL DNLS 
equation within the framework of the bilinear formalism. In particular, we present two 
different expressions of the bright iV-soliton solution in a simple determinantal form. The 
proof of the solution can be obtained by a straightforward calculation where one relies 
only on a few basic formulas for determinants. This paper is organized as follows. In 
section 2, we develop a direct method of solution. Specifically, we first transform equation 
(1.1) to a system of bilinear equations by introducing an appropriate dependent variable 
transformation. Notably, we show that one of the bilinear equations can be replaced by 
a trilinear equation which is easy to prove when compared with other bilinear equations. 
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In section 3, we present the bright iV-soliton solution of the bilinear equations. It has a 
simple structure expressed in terms of certain determinants. Subsequently, we perform 
the proof of the bright iV-soliton solution using an elementary theory of determinants in 
which Jacobi's identity will play a central role. Last, we provide an alternative expression 
of the bright iV-soliton solution. The equivalence of both expressions is discussed by 
referring to the properties of the Cauchy matrix. We emphasize that the bright 7V-soliton 
solution obtained here yields the bright iV-soliton solution of the derivative NLS equation 
by replacing simply the time dependence of the solution. This connection has also been 
found in the context of the 1ST [4]. In section 4, we investigate the properties of the 1- 
and 2-soliton solutions and then analyze the asymptotic behavior of the bright iV-soliton 
solution in the limit of infinite time. As a result, we obtain the explicit formula for the 
phase shift in terms of the amplitude parameters of solitons. Section 5 is devoted to 
concluding remarks. 

2. Exact method of solution 

In this section, we develop a direct method of solution for constructing bright soliton 
solutions of equation (1.1) under the boundary condition u — > as \x\ — > oo. Specifically, 
we show that equation (1.1) can be transformed to a system of bilinear equations. We also 
demonstrate that one of the bilinear equations can be replaced by a trilinear equation. 

2.1. Bilinearization 

The bilinearization of equation (1.1) is established by the following proposition: 
Proposition 2.1. By means of the dependent variable transformation 



equation (1.1) can be decoupled into the following system of bilinear equations for the tau 
functions f and g 




(2.1) 



D x D t g ■ f = gf, 



(2.2) 



DJ ■ f* = igg*, 



(2.3) 



D x D t f-f* = iD x g-g*. 



(2.4) 
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Here, f = f(x,t) and g = g(x,t) are complex-valued functions of x andt, and the asterisk 
appended to f denotes complex conjugate and the bilinear operators D x and D t are defined 
by 

( d d \ m f d d\ n 



where m and n are nonnegative integers. 



Proof. Substituting (2.1) into (1.1) and rewriting the resultant equation in terms of the 
bilinear operators, equation (1.1) can be rewritten as 



jz(D x D t g ■ f - gf) - j§- 



2 (D x D t g -f-gf)- -§-(f*D x D t f ■ f + 2ig*D x f -g)=0. (2.6) 



Inserting the identity 

f*D x D t f ■ f = fD x DJ ■ f* - 2f x D t f ■ f* + /(A/ • f%, (2.7) 

which can be verified by direct calculation into the second term on the left-hand side of 
(2.6), one modifies it in the form 

j 2 {D x D t g -f-gf) 
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■[f(D x D t f ■ f* - iD x g ■ g*) - 2f x (D t f ■ f - igg*) + f(D t f ■ f* - igg%} = 0. (2.8) 
By virtue of equations (2.2)-(2.4), the left-hand side of (2.8) vanishes identically. □ 

It follows from (2.1) and (2.3) that 



12 , 9 ,_ /* 



The above formula gives the modulus of u in terms of the tau function /. 
2.2. Remark 

Proposition 2.2. The trilinear equation for f and g, 

f*(g x tf - 9th - gf) = fligJ - gU). (2.10) 

is a consequence of the bilinear equations (2.2)-(2.4)- 



4 



Proof. Using the definition of the bilinear operators, one can show that 

f{9xtf ~ 9tfx ~ 9f) - ft(9xf ~ 9 fx) 

= f*(D x D t g. f - g f)- g -{(D t f ■ f -igg% + (D x D t f • f -iD x g- g*)} + g x (D t f ■ f -igg*). 

(2.11) 

The right-hand side of (2.11) becomes zero by (2.2)-(2.4), which yields (2.10). □ 

In view of proposition 2.2, we can use equation (2.10) in place of (2.2), for example in 
performing the proof of the bright iV-soliton solution. We will see later in section 3.6 that 
the proof of (2.10) is simplified considerably compared with that of (2.2), even if (2.10) 
in not bilinear but trilinear in / and g. 

3. Bright iV-soliton solution and its proof 

In this section, we show that the tau functions / and g representing the bright 7V-soliton 
solution admit the compact determinantal expressions. This statement is proved by an 
elementary calculation using the basic formulas for determinants. Then, we provide an 
alternative form of the bright iV-soliton solution. Last, we demonstrate that the tau 
functions / and g satisfy a trilinear equation. The implication of this result will be 
discussed briefly in conjunction with the bright iV-soliton solution of the DNLS equation. 

3.1. Bright N-soliton solution 

The main result in this paper is given by the following theorem. 

Theorem 3.1. The bright N-soliton solution of the system of bilinear equations (2.2)- 
(2.4) is expressed by the following determinants 

f=\D\, g = 

Here, D is an N x N matrix and z, z t and 1 are N -component row vectors defined below 
and the symbol T denotes the transpose: 

z ■ z* in* / \ \ 

D = {d jk )i<^ k < N) d jk = 3 k Zj = exp pjx H 1 + Qo , (3.2a) 

Pj + Pk V Pj J 



D 



(3.1) 
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z = (z!,z 2 , ...,z N ), z t = (z 1 /p 1 ,z 2 /p 2 , z N /p N ), 1 = (1, 1, 1), (3.26) 
where pj and Q (j = 1, 2, N) are arbitrary complex parameters. 

The bright iV-soliton solution is parameterized 2N complex parameters pj and Qo (j = 
1,2,...,N). The parameters pj determine the amplitude and velocity of the solitons 
whereas the parameters (j determine the phase of the solitons. The conditions pj +p* k ^ 
for all j and k and pj ^ pt for j ^ k may be imposed on the parameters to assure the 
regularity of the solution. We point out that the iV-soliton formula given by Theorem 
3.1 can be shown to take the same form as that obtained in [4] by changing the soliton 
parameters pj and Q (j = 1, 2, N). 

3.2. Notation and basic formulas for determinants 

First, we define the definition of matrices associated with the bright iV-soliton solution 
and then provide some basic formulas for determinants. The following bordered matrices 
appear frequently in our analysis: 

(D c T d T \ 
a j , (3.3) 

where a, b, c and d are iV component row vectors. Let Djk be the cofactor of the element 
djk- The following formulas are well known in the theory of determinants: 

j,k=i 

D a T A 

b = \D\z - 22 D jk ajb k , (3.5) 
j,k=i 

|L>(a,b;c,d)||L>| = |D(a;c)||D(b;d)| - |D(a;d)||D(b;c)|. (3.6) 

Formula (3.4) is the differentiation rule of the determinant and (3.5) is the expansion 
formula for a bordered determinant with respect to the last row and last column. Formula 
(3.6) is Jacobi's identity. 

3. 3. Differentiation rules and related formulas 
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In terms of the notation (3.3), the tau functions / and g can be written as / = \D\ and 
g = |-D(1; z t )|, respectively. The differentiation rules of the tau functions with respect to 
t and x are given by the following formulas: 

Lemma 3.1. 

f t = -\Dte;zt)\, (3.7) 

f x = -\D(z*;z)\, (3.8) 

f xt = -\D(z* t ;z)\ - \D(z*;z t )\ + |D(z*, z*; z, z t )|, (3.9) 

<7 t = |£>(l;z„)|, (3.10) 

g x = \D(l;z)\-\D(l,z*;z t ,z)\, (3.11) 

^ = p(l;z t )|-p(l,z*;z tt ,z)|. (3.12) 

Proof. We prove (3.7). Applying formula (3.4) to / given by (3.1) with (3.2a), one 
obtains 



N 

ft = D * 



Z j Z k 



hk=l P 3 Pk 
N 

j,k=l 

where in passing to the second line, use has been made of the relations Zj/pj = Zj tt , z* k p* k = 
z* k ^ t . Referring to formula (3.5) with z = and taking into account the notation (3.3), 
the above expression becomes the right-hand side of (3.7). Formulas (3.8)-(3.12) can be 
proved in the same way if one uses (3.4), (3.5) and the relation z xt = z as well as some 
basic properties of determinants. □ 

The complex conjugate expressions of the tau functions / and g and their derivatives are 
expressed as follows: 

Lemma 3.2. 

f* = \D\, D = D* T , (3.13) 
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/ t * = -|L>(z t *;z t )|, (3.14) 

r x =-\D(z*;z)\, (3.15) 

f* xt = -\D(z*;z t )\ - \D(z* t ;z)\ + \D(z* : z t *; z, z t )|, (3.16) 

<f = |5(zT;l)|, (3.17) 

i/: = |5(z*;l)|-|5(zr,z*;l,z)|, (3.18) 

where D = (dj k )i<j t k<N is an N x N matrix with elements dj k = dj k + i. 

Proof. It follows from (3.2a) that d* k = d k j + i or equivalently, D* = D T , which proves 

(3.13) since /* = \D*\ = \D T \ = \D\. Formulas (3.14)-(3.18) can be proved in the same 

way. □ 
The formulas below are used to reduce the proof of the bright iV-soliton solution to 
Jacobi's identity. 
Lemma 3.3. 

|£>(l;z)| = |D(l;z)|, (3.19) 

|D(l;z t )| = |L>(l;z t )| = (-l^^p*; z)|, (3.20) 

\D(z*;z)\ = (-l) N ~ 1 c\D(z* t ;z x )\, (3.21) 

\D(z*;z t )\ = (-l) N - 1 c\D(z* t ;z)\, (3.22) 

\D(z*;z tt )\ = (-lf- 1 c\D(z* t ;z t )\, (3.23) 

\D(z*, 1; z, p)| = (-l) N c\D(zl p*; z x , 1)|, (3.24) 

|L>(l,z*;z,z t )| = \D(l,z* t ;z,z t )\, (3.25) 

where p = (1/pi, I/P2, 1/pat) z x = (piZi,p 2 Z2, Pn%n) (ire N component row 
vectors and c = Ylf=i(P*j / Pj) ■ 

Proof. We prove (3.21) only. The other formulas can be proved following the same 
procedure as that described below. First, using the key relation which follows from (3.2a) 
and the definition of djk, 

d jk - Z &- = - P l{d jk + i) = - P M jk , (3.26) 
Pj Pj Pj 
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the left-hand side of (3.21) can be modified into the form 

|L>(z*;z)| = 



z* 







(3.27) 



after multiplying the (N + l)th row by —Zj/pj and then adding the resultant expression 
to the jih row for j = 1,2, ...,7V. Formula (3.21) follows from (3.27) if one extracts the 
factors 1/pj and — p* k from the jth row and kth column, respectively for j, k = 1, 2, N 
and the factor —1 from the (N + l)th row. □ 

3.4- Proof of the bright N-soliton solution 
3.4.1. Proof of (2.2). Let 

Pi = D x D t g ■ f - gf = g xt f - g x f t - g t f x + gf xt - gf. (3.28) 

Substituting formulas (3.7)-(3.12) into this expression and applying Jacobi's identity to a 
term \D(1, z*; z tt , z)\\D\, Pi reduces to 



Pi — Pll + Pl2 + Pl3: 



with 



P 11 = \D(l;z)\\D(z* t ;z t ) 
P 12 = \D(l;z t )\\D(z*,z* t ;z,z t ) 
P 13 = \D(l;z)\\D(z*;z tt ) 
Referring to Jacobi's identity, Pn becomes 



-|£>(l;z t )||L>(z t *;z)|, 

-\D(z* t ;z t )\\D(l,z*;z t ,z)\, 

-\D(l;z t )\\D(z*;z t )\. 



Consider the identity 







Pu 


= \D\\D{\, 


* 




\D 


(i; 


*)l 


l^(i; z*)l 


\D 


(i;*)l 


\D( 




*)l 


I^K;z t )| 


\D( 


z t *;z)| 


\D( 




*)l 


\D(z*;z t )\ 


\D( 





(3.29a) 

(3.296) 
(3.29c) 
(3.29d) 

(3.30) 



0, 



(3.31) 
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which follows immediately since the first two columns of the determinant coincide. Ex- 
panding the above determinant with respect to the first column and using Jacobi's identity, 
one finds that 

P 12 = p(z*;z t )||D(l,z t *;z,z t )|. (3.32) 
Gathering up three terms (3.29d), (3.30) and (3.32), Pi simplifies to 

(3.33) 

It follows from formula (3.5) that 



\D\ + \D(z*;z t ) 



D zj 
z* 1 



-l) N c\D\, 



(3.34) 



where in passing to the last line, the procedure used for deriving lemma 3.3 has been 
applied. If one substitutes this result together with (3.19), (3.20), (3.23) and (3.25) into 
Pi, one can recast Pi into the form 

Pi = (-l) N - l c{-\D\\D(l, z*; z, z t )\ + |P(1; z)||P(z*; z t )\ - \D(1; z t )\\D(z* t ; z)|}, (3.35) 



which becomes zero by virtue of Jacobi's identity. 
3.4.2. Proof of (2.3). Let 



P2 = D t f ■ f* - igg* = ftf* - fft - igg*- 



□ 



(3.36) 



Substituting formulas (3.7), (3.13), (3.14) and (3.17) into this expression, P 2 becomes 



P 2 = -|P>(z t *;z t )||P>| + |P>||P>(z t *;z t )| - i\D(z*; l)\\D(l;z t ) 



(3.37) 



The formulas below can be derived by using the definition of the matrix D from (3.13): 



\D\ = \D\-i\D(l;l)\, 

|P(z-l)| = |P(z-l)|, 
\D(z* t ;z t )\ = |P>(z t *;z t )|-i|P>(z t *,l;z t ,l) 



(3.38) 

(3.39) 
(3.40) 
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After introducing (3.38)-(3.40) into (3.37), P 2 reduces to 

P 2 = i{\D(z* t ; z t )\\D(l; 1)| - \D\\D(z* t , 1; z t , 1)| - |£>(z t *; 1)||£>(1; z t )|}. (3.41) 

One sees that P 2 becomes zero by Jacobi's identity. 
3.4-3. Proof of (2.4). 

Instead of proving (2.4) directly, we differentiate (2.3) by x and add the resultant expres- 
sion to (2.4) and then prove the equation P 3 = 0, where 

P 3 = htf* - fj; - ig x g*. (3.42) 

Substituting (3.8), (3.9), (3.11), (3.13), (3.14) and (3.17) into (3.42), P 3 becomes 

P 3 = {-\D(z* t ;z)\ - \D(z*;z t )\ + \D(z*, z* t ; z, z t )\)}\D\ - \D(z*;z)\\D(z* t ;z t )\ 

-i{\D(l;z)\-\D(l,z*;z t ,z)\}\D(z* t ;l)\. (3.43) 
Using the cofactor expansion of the determinant \D\, 

N N 

^d jk D jk = \D\, ^d jk D jk = \D\, (3.44) 

k=l j=l 

one obtains 

V (- + \) d jk D jk = V (- + \) \D\. (3.45) 

The definition of dj k from (3.2a) gives ( — + 4r ) d jk = z^ t z* kt — ipj, which, after substi- 
tuting into (3.45) and using formula (3.5), yields the relation 

N ( i n 

-\D{z\- z t )\ + i|D(l; p)| = - + "J )\D\. (3.46) 

j=i \Pi PjJ 

Invoking the relation D* = D T , the complex conjugate of this expression can be written 
in the form 

N 



D(z* t ; z t )\ - i\D(p*; 1)| = £ (I + I) \D\. (3.47) 

j=1 \Pj PjJ 
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The following formula can be derived if one differentiates (3.46) by x and notes the relation 

Zxt = Z 

-\D(z* t ;z)\ - \D(z*;z t )\ + \D(z*, z* t ; z, z t )\ - i\D(z*, 1; z, p)| 



N 



Multiply (3.48) by \D\ and use (3.47) to obtain 

{-\D(z* t ;z)\ - \D(z*;z t )\ + \D(z*, z* t ; z, z t )\}\D\ 

= i|^l;z,p)|p| + {|D(z:;z t )|+i|5(r;l)|}P(z*;z)|, 
which, substituted into P 3 , gives 

A = i 



(3.48) 



(3.49) 



|D(z%l;z,p)||D| + |D(p*;l)||D(z-z)|-{|D(l;z)|-|D(l,z-z,,z)|}|D(z-i; 

(3.50) 

Differentiation of (3.20) with respect to x yields the relation 



|D(l;z)| - \D(l,z*;z t ,z)\ = (-lf- l c\D(p*;z x )\. 



(3.51) 



If one substitutes (3.21), (3.24) and (3.51) into (3.50), P 3 simplifies to 

P 3 = i(-l)^{|£(z^;z x ,l)||£H (3-52) 

It turns out that P 3 becomes zero by Jacobi's identity. □ 
3. 5. An alternative expression of the bright N-soliton solution 

Here, we present an alternative expression of the bright iV-soliton solution of equation 
(1.1). It is expressed by the following theorem: 

Theorem 3.2. The tau functions f and g' given below satisfy the system of bilinear 
equations (2. 2)- (2.4): 



r 



A I 
-I B 



9 = 



A I yf 
-I B T 
10 



(3.53) 
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Here, A, B and I are N x N matrices and y and y t are N component row vectors defined 
respectively by 

A = {a jk )i< jik < N) a jk = - ' \ , y-j = exp I qjX H t + rj j0 , (3.54a) 

5j + Ik \ 5j J 

B = {bjk)i<j,k<N, b jk = — , (3.546) 

Qj + Qk 

I = (Sjk)i<j,k<N : N x N unit matrix, (3.54c) 

y = (yi,V2, -,vn), y* = (2/1/51,2/2/52, -,W5iv)- (3.54d) 

where qj and rjj (j = 1,2, ...,N) are arbitrary complex parameters. 

The proof of theorem 3.2 parallels theorem 3.1 and hence omitted. Instead, we provide 
an alternative proof. To this end, we first establish the following proposition: 

Proposition 3.1. Under the parameterization qj = —p*j,i]jo = — C* +ln Cj (j = 1, 2, N), 

the tau functions f,g, f and g' satisfy the relations 

f = (-1)^1/, g' = -c'\A\g, (3.55) 

where 

and the conditions qj + ql ^ for all j and k and qj ^ q k for j ^ k are imposed on the 
parameters. 

Proof. By means of the operation of matrix multiplication, one can show that 

f = \I + AB\ = \A\\A- 1 + B\, g'= 1 y J = \A\ A ~' + B . (3.57) 

The inverse of the Cauchy matrix A exists due to the conditions imposed on the param- 
eters. It reads [7] 

A"' = ( ..!. ) ■ (3^8) 



q* + q k y*y k/ ^ k < N 
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The specified parametrization also leads to the relation yj = Cj/z*, which, substituted 
into A^ 1 and B, gives 

Pj + Pk) l<j,k<N 



A- l = 



(3.59) 



J>3 + Pk) l<j,k<N 

It follows from (3.2a) and (3.59) that A^ 1 + B — —D. If one introduces this relation into 
/' from (3.57), one finds the first relation of (3.55). To proceed, note that the j'th element 
of the column vector A~ 1 yf is given by 

[1 

^ )jmy m ,t ^ n _ri 

m=l 



N 



(^- i yn i =E^" 1 w™.*=SE 



i nii(ft»+«r) 



y'j t^i 9"*Ilj=i x (9m 



(3.60) 



jv-i rr^ 



The sum with respect to m turns out to be {—l) N ~ 1 Yli = i{qi/qi)/q] by using Euler's 
formula. Hence, (3.60) reduces to (— l) N c'zj tt upon rewriting the resultant expression in 
terms of Zj. Thus, g' from (3.57) becomes 



-c'\A\g, 



9 = (-lfc\A\ | ° ^ 

giving rise to the second relation of (3.55). □ 

Proposition 3.2. The tau functions f and g' from (3.55) satisfy the system of bilinear 
equations (2.2)-(2.4) if f and g satisfy the same system of equations. 

Proof. Let P[, P' 2 and P' z be 

P[ = D x D t g' ■ f - g'f = gjf - g'J' t - g'J' x + g'f xt ~ f , (3-61) 

PL = D t f ■ f* - igfgT = fX ~ f'f't - tig'*, (3.62) 

^3 = fLX - fXt - WJ*, (3.63) 
respectively. First, note that the determinant \A\ from (3.54a) can be modified in the 
form 

N N 

\A\=Uy s ]lyl 



3=1 k=l 
N 



Qj + % 



= exp 



jo + Vjo 



1 



Qj + 9* 



(3.64) 
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It immediately follows from this expression that 

D x D t \A\ ■ \A\ = 0. (3.65) 
Substitute (3.55) into (3.61)-(3.63) and use the relation |c'| = 1 to obtain 

P[ = (-lf-'c'ilAWD^g ■f-gf) + (D x D t \A\ ■ \A\)gf}, (3.66) 

P^ = \A\ 2 (D t f-f*-igg*), (3.67) 

P> = \A\\A\ X (DJ ■ f* ~ igg*) + |^| 2 (/, t r - fJt* ~ ^9*) + \{D x D t \A\ ■ \A\)ff. (3.68) 

Above three expressions vanish identically by virtue of (2.3), (2.4) P 3 = with P 3 given 
by (3.42), and (3.65). □ 

Propositions 3.1 and 3.2 lead to an alternative expression of the bright iV-soliton 
solution of equation (1.1) in terms of the tau functions /' and g': 

u=(-lf- 1 c' 9 ? , H 2 = i|ln^. (3.69) 

Note that if u satisfies equation (1.1) then cu with c being a complex constant satisfies 
the equation as well if |c| = 1. Hence, the factor (— l^^V in the above expression for u 
is irrelevant and it may be replaced simply by 1. 

3.6. Remark 

As already mentioned in section 2.2, one of the bilinear equations (2.2)-(2.4) can be 
replaced by the trilinear equation (2.10). Here, we show that the tau functions (3.1) for 
the bright iV-soliton solution satisfy equation (2.10). The proof is quite simple. Indeed, 
referring to (3.7), (3.8), (3.10)-(3.14) and (3.23) and using Jacobi's identity, we can derive 
the following relations 

g x tf - gtf* -gf = (-1)^1^(1; z*)l = (-1)^(1; z)\f t *, (3.70) 

9j-9f x = (-l)"c|£>(l;z)||5| = (-lfc\D(l;z)\r. (3.71) 

Upon substituting (3.70) and (3.71) into (2.10), one can confirm that equation (2.10 ) 
holds identically. Note that expressions (3.70) and (3.71) have a common factor \D(1; z)|. 
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Let q = (g/f) x - K then follows from (3.71) that 

q=( _ lfc m^r (3 . 72) 

The form of q coincides perfectly with the iV-soliton solution of the DNLS equation 

iq t + q xx + 2i{\q\ 2 q) x = U. (3.73) 

Actually, if we replace the time dependence of Zj in (3.2a) as ip|i instead of (l/pj)t, 
then the iV-soliton solution (3.72) with the tau function / from (3.1) satisfies the DNLS 
equation. See, for example [8-10]. For completeness, we reproduce the 1-soliton solution 
qi as well as the square of its modulus |gi| 2 . To this end, we put p± — a± + ibi and 
Cio — $io + iXio- Then, the expression (3.72) with N — 1 yields 

2a 1 (a 1 -io 1 ) e^je^-h + ia,) 

ft = ^7 T~w, 1 ~~Y2 — ' 3.74a 

ai + i6i (e 2til — bi — iai) 2 

Ox = ai(x - 2bit) + 0io, Xi = hx + {a{ - b\)t + X io, (3.746) 

9n 2 1 1 
2 = Za l i e 2<5i = 1 ^ 

xAiTfef cosh 2(9 1 + 5 1 )- ' v^Tfe!' 

Let Ai and c\ be the amplitude and velocity of |gi|, respectively. We then find from (3.75) 
that 



Al = V V c i + 4a i + Cl ' ci = 26i. (3.76) 
Note that if foi > 0, the soliton propagates to the right whereas if bi < 0, it propagates to 
the left. This fact implies that both the overtaking and head-on collisions are possible in 
the interaction process of two solitons. It is interesting that in the case of bi > 0(&i < 0), 
the amplitude of the soliton is an increasing (a decreasing) function of the velocity. 

The result (3.72) also follows from the view point of the 1ST. Indeed, if we identify q 
with u x , then the x-part of the Lax pair for the FL DNLS equation coincides with that 
of the DNLS equation [2]. It turns out that the bright iV-soliton solution of the DNLS 
equation can be derived from that of the FL DNLS equation via the relation q = u x and 
replacing the time dependence of the exponential functions Zj, as demonstrated here. 



16 



4. Properties of the bright iV-soliton solution 



In this section, we investigate the properties of the bright iV-soliton solution particularly 
focusing on the 1- and 2-soliton solutions and then address the asymptotic behavior of the 
iV-soliton solution for large time. To this end, we first parametrize the complex constants 
Pj and Qo by the real constants a,, bj, 9j and Xjo as 

Pj = a j + ib ji Go = Ojo + iXjo, aj > 0, j = 1, 2, N, (4.1) 

and introduce the new independent variables 9j and Xj m accordance with the relations 

1 

0j = aj(x + Cjt) + j0 , Cj= 2 2 , j = l,2,...,N, (4.2a) 

a j + °j 

Xj = bj{x - Cj t) + Xj o, J = 1,2,..., N. (4.26) 
In terms of these variables, Zj is expressed as 

Zj = e°i +i *, j = 1,2, ...,N. (4.2c) 

^.i. Bright 1-soliton solution 

The tau functions f — fi and g = gi corresponding to the bright 1-soliton solution follows 
from (3.1) and (3.2) with N = 1. They read 

/. = 91 = - Z ~ (4-3) 

Pi + pj Pi 

In terms of the new parameters defined by (4.1) and (4.2), the 1-soliton solution u\ = 
u i(9i,Xi) takes the form of an envelope soliton 

Ul = ^_1±A *i = e9l+iXl (4.4) 

Pi z\z\ — \p\ a\ + ifoi e 2 " 1 — b\ — \a\ 

The square of the modulus of U\ can be written as 

| Ul I 2 = ^1 e 251 = 1 U 5a) 

(a 2 + 6 2 )f cosh 2(0 1 + 6 1 )--j!^' yffi+W 
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showing that \u±\ has the amplitude A-y given by 



(4.56) 



V^T~bJ 

It represents a localized pulse (or a bright soliton) moving to the left at a constant velocity 
ci = \/{a\ + b\). The amplitude- velocity relation follows immediately from (4.5b) to give 



A 1 = + 



(4.6) 



See Figure 1. The characteristic of the solution depends crucially on the signature of the 
parameter b±, as we shall see now. 

2 
1.5 

< 1 
0.5 




0.2 0.4 0.6 

Ci 



0.8 1 



Figure 1. Amplitude- velocity relation for bi — 1 (solid line) and b± — — 1 (broken line). 

For a fixed value of bi, Ai is an increasing function of c\ if b\ > 0. If, on the other 
hand, b\ < 0, then Ai becomes an increasing function of c\ in the interval < c\ < -^j 
and a decreasing function of C\ in the interval ^ < C\ < ^. In the case of b\ > 0, the 
amplitude A\ remains finite in the limit of a\ — > (i.e., infinite width) as opposed to the 
behavior of the usual soliton for which the amplitude tends to zero in the limit of infinite 
width. It then turns out that the bright soliton u\ from (4.4) reduces to an algebraic 
soliton of the form 



lcf exp 



Ul = 



x + c\t + x — i^ 
18 



C l = 12' 



1 



(4.7a) 



and |wi| 2 from (4.5) becomes 



M = 7 — (4-76) 

where we have put #io = airr — 5i. The similar peculiar structure of the solution has 
been found for the bright soliton solution of the DNLS equation [11]. 

In the case of bi < 0, we can observe that when the velocity lies in the interval 
< ci < ^j, the solution exhibits an interesting feature. Specifically, a large soliton 
propagates slower than a small soliton, as will be illustrated graphically in the interaction 
process of two bright solitons (see figure 3). The amplitude- velocity relation described 
above is in striking contrast to that of the 1-soliton solution of the DNLS equation. 
Actually, we see from (3.76) that A\ becomes a decresing function of |ci| when b\ < 0. 

Last, it is instructive to compare the solution (4.4) and (4.5) with that derived by the 
1ST. We define the new parameters A x and 71 by the relations a\ = A 2 sin 71, b x = 
—A 2 cos 71, < 71 < 7T, and then replace 9 W by 9 W + InAi (or equivalently, Q\ by 
Q\ + In A x ). We find that expressions (4.4) and (4.5) recast respectively to 

1 2ie ei+ixi sin 71 .„ , 

Ul = "AT c^+m + 1 ' (48a) 

2 1 2 sin 2 71 
H = Af cosh 2^ + cos 71 ' ( } 

which is just the 1-soliton solution presented in [3]. Note in these expressions that A 1 = 
2 sin^/Ai and c\ = 1/A 2 . The algebraic soliton is generated from (4.8) if one puts 
71 = 7r — e and takes the limit e — > +0 [3] . 

4-2. Bright 2-soliton solution 

The bright 2-soliton solution exhibits a variety of interesting features. Here, we investi- 
gate the interaction process of bright solitons focusing on the asymptotic behavior of the 
solution. As a result, we obtain the formulas for the phase shifts for each soliton. The 
corresponding tau functions f 2 and g 2 are given by (3.1) and (3.2) with N = 2. Explicitly, 

, = plpli.pi -p 2 )(pi -p* 2 ) z z * 

{P1+P\){PI+PI){P2+P\){P2+PI) {Pl+P\){P2+P* 2 ) ZlZl 
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* i 



ipl 



(Pi + P*l)(P2 + P* 2 ) (Pi + P* 2 )(P2 + P*l) 

(pi -P2)(pl -p* 2 ) 



Z\Z 2 



(pi + p\) (pi + Pi) (p 2 + Pi) (P2 + P2) 



#2 = 



ip 2 (p* -p*) 



-21 + 



ipi(Pi -P2) 



-22 



(pi+P2)(P2+Pl) 2 1 

Z±Z 1 Z 2 Z 2 , 

p*l(pl -P2) 



(4.9a) 



Pl(P2 +P*l)(P2 + P 2 ) Pl{V\ +P*l)(Pl +P* 2 ) PmiPl +P\){P2 +Pi) 

P 2 (Pl~P2) * 
—Z\Z 2 Z 2 - 



Z\Z X z 2 



(4.96) 



P\P2{P\ +P* 2 )(P2 +P* 2 ) 

Now, we order the magnitude of the velocity of each soliton in the (x, t) coordinate 
system as C\ > c 2 . We first take the limit t — > —00 with 9 1 being fixed. Since in this limit 
\z\ I =finite and \z 2 \ — > 00, the leading-order asymptotics of f 2 and g 2 are found to be as 

(Pi -P2W1 ~P* 2 ) 



/a 



^2^2 

(pi +Pt)(p 2 + P2) KPl+PsDfo+Pl) 
P2(Pl -P2 



- ip x 



(4.10a) 



92 



-Z\Z 2 Z 2 . 



PlP2(Pl + P* 2 )(P2 + P* 2 ) 

The asymptotic form of the 2-soliton solution u 2 evaluated by (4.10) is expressed as 



(4.106) 



u 2 



Pi +P\ 



Pi z[z[* -\pi 



z 1 = Z\ exp 



-In 



P2(Pl +P2) 



(4.11) 



.P* 2 (P1 -P2) 

Thus, the asymptotic of u 2 takes the same profile as that of the 1-soliton solution (4.4) 
except the phase shifts. Specifically, 



U2~«i(0i + A0j \ X i + ^x\ 



A9\ 



(-) 



-In 



Pi+Pl 



(-) 



Pi +P 2 P2 . 

- arg arg — + n. 

Pi - P2 P 2 



(4.12a) 
(4.126) 



Pi -P2 

Next, we take the limit t — > +00 with 6\ being fixed. In this limit, \z\\ =finite and 
\z 2 \ — > 0. The expressions corresponding to (4.10) and (4.11) are given by 



ip* 2 {pi-p 2 ){p\-p* 2 ) 



{pi+P* 2 ){p 2 +p\) 
(pi+P* 1 )(pi+P 2 )(p 2 +P* 1 )(p 2 +P 2 ) { (pi -p 2 )(Pi -P* 2 ) 

m(p*!-p* 2 ) 



z x z x - ip l 
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u 2 



Pl+Pl 



II ff* 

Pi z'[z\ -lpi 



Pl{P2 + P\){P2 + P* 2 ) 

z'l = z\ exp 
20 



Zl, 



ln(^±4 
{p* 2 {pi-p 2 ) 



(4.13a) 

(4.136) 
(4.14) 



It follows from (4.14) that 



M2 ~ Ml (^ + A^ +) , Xl + A x i 



A^ +) = In 



Pl + P*2 



Pl ~P2 



a (+) P1+P2 , P2 , 

Axi = arg + arg — + vr. 

Pi - Pi P2 



(4.15a) 
(4.156) 



In view of the fact that solitons propagate to the left and the trajectory of the center 
position of the soliton is described by the equation Q\ + A9^ = 0, we can define the 



phase shifts Axj = (AO^ — A9j~')/aj and Axj = ^xf ~ ^-Xj~' which represent the 
total changes of the center position and the phase of the jth soliton (j = 1, 2), respectively. 
Using (4.12) and (4.15), we then obtain 



,(+) 



Ax 1 = — In 
ai 



P1+P2 



P1-P2 



Axi = 2 arg 



P1+P2 , P2 

h 2 arg— . 

Pi ~ P2 P 2 



(4.16) 



We can perform the similar asymptotic analysis while keeping 62 fixed. Hence, we quote 
the final result. As t — > —00, the expression corresponding to (4.12) reads 

u 2 ~ Ul (9 2 + A9i~\ X 2 + A X ( 2 ] ), 



A6 { f ] = In 



P2 + P\ 



P2+Pl , Pl , 

arg h arg — + ir, 

P2 - Pi p\ 



P2-P1 

whereas as t — > +00, the expression corresponding to (4.15) reads 



u 2 ~ ui(Q 2 + A#2 +) , X2 



Ad { + ] 



-In 



P2+P1 



P2 ~Pl 



Ax 



(+) 



— arg 



-ax { 2 +) : 
P2+PI 

P2 ~Pl 



Pl , 

arg — + 7r. 

Pi 



Thus, the total phase shifts are given by the formulas 



Ax 2 = In 

a 2 



P2+PI 



P2 ~Pi 



AX2 = -2 arg 



P2+PI 

P2 ~Pl 



2 arg— . 
Pi 



(4.17a) 
(4.176) 

(4.18a) 
(4.186) 

(4.19) 



Figure 2 shows the time evolution of a bright 2-soliton solution U — \u\ with positive 
bj (j = 1,2). See the solid line in fugure 1 which plots the amplitude- velocity relation for 
bi = 1. In the present example, A 1 = 1.99, A 2 = 1.44, ci = 0.99, c 2 = 0.20 and the phase 
shifts are given by Axi = 2.00, Ax 2 = —0.10. The feature of the interaction process of 
solitons is similar to that of the usual solitons, namely the large soliton propagates faster 
than the small soliton and the former one suffers a positive phase shift whereas the latter 
one suffers a negative phase shift. 
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Figure 2. Time evolution of a bright 2-soliton solution with the parameters a\ = 0.1, 02 = 
2.0,6i = 1.0, b 2 = 1.0. 

Figure 3 depicts the similar plot for negative bj (j = 1, 2) where A\ = 0.43, A 2 = 
0.69, ci = 0.80, c 2 = 0.41 and Axi = 3.55, Ax 2 = -1.48. 




Figure 3. Time evolution of a bright 2-soliton solution with the parameters a\ = 0.5, a 2 = 
1.2,6! = -1.0,6a = -1.0. 



22 



In this case, the velocities Cj(j = 1, 2) are chosen to satisfy the inequality < Cj < ^ 



Consequently, the velocity of each soliton becomes a decreasing function of the amplitude, 
as indicated by the broken line in figure 1. Unlike the previous example, there appears a 
peculiar feature in the interaction process of solitons. Indeed, as evidenced from figure 3, 
the large soliton is seen to propagagate slower than the small soliton and suffers a negative 
phase shift. On the other hand, the small soliton suffers a positive phase shift. 

4-3. Bright N- soliton solution 

The asymptotic analysis for the bright iV-soliton solution can be performed straightfor- 
wardly following the 2-soliton case. Hence, without entering into the detailed calculation, 
we describe the main results only. To proceed, let us order the magnitude of the velocity 
of each soliton as c\ > C2 > ... > cn- We take the limit t — > — oo with 9 n being finite. Since 
in this limit \zj\ — > for j < n and \zj\ — > oo for n < j, we find that the leading-order 
asymptotic of the tau function / = f N from (3.1) and (3.2a) can be written in the form 

n-l N 

f N ~d{l,2,...,n)d{n+l,n + 2,...,N) JJ(-ip fe ) (z jZ *) (d n z z z* n - ip* n ), (4.20a) 

k=l j=n+l 

where 



d(m, m+1, n) 



J>3 + PkJ m<j,k<n 

is the Cauchy determinant and 



Um<j<k<n(Pj - Pk)(Pj - Pi) 
Y[m<j,k<n(Pj + P*k) 



1 < m < n < N, 
(4.206) 



T-rN (Pn-Pj)(Pn~Pj) 

_ d(l,2,...,w- l)d(n,n + l,...,N) _ lij=n+i ( Pn +P*)( P * n + Pj ) ^ 
n ~ d(l,2,...,n)d(n+ l,n + 2,...,N) ~ ri^ (p^p^p^p*) ^ " c > 

llj=l (Pn+P*)(pl+Pj) 

Similarly, has the asymptotic form 

(-1) N 

g N ~ ^ J —d(l,2,...,n-l)d(n + l,n + 2,...,N) 

Pn 

*f\(-Vk)T\ Pi( f n ~ P *l TT PliPn ~ Pk ] IT (zjz*)z n . (4.21) 
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Taking into account the relation 



d(l,2,...,n) 



n-l 



n 



{Pn-PjWn-p*) 



d(l,2,...,n-l), 



Pn + Pn | = i (Pn + Pj)(Pn + Pj) 

we obtain from (4.20) and (4.21) the asymptotic form of the iV-soliton solution 



(4.22) 



U N 



Pn fJl P* iPn ~ Pj ) k *} +1 Pk iPn + P* k ) d n Z n Z* - ip* 



This expression can be rewritten in terms of the 1-soliton solution as 

u N ~ Ul (9 n + A6(-\ Xn + A X i _) ), 

with 



n-l 

A xi _) = \ arg 
j'=i 



n-l 

Ae ) = E in 

i=i 

Pn+P* 



Pn + P* 



Pn - Pj 



+ arg — 



Pn ~ Pj 



Pj 



N 

- E 

j=n+l 
N 

E 

j=n+l 



Pn+P*i 



arg 



Pn - Pj 

Pn + P*j 
Pn - Pj 



Pi 

+ ar § ( "T 

\P* 



As t — )> +oo, on the other hand, the asymptotic form of un is found to be as 
, N Pn +P n TT + Pj) TT ~ 2 ™ 



MAT ~ (-1) 



n 



n 



(4.23) 

(4.24a) 
(4.246) 

(4.24c) 
(4.25) 



Pn M +1 P*j (Pn - Pj) ^ Pk(Pn + Pfc) d n ^n< ~ 1P„ ' 

which leads to the expression 

u N ~ «i(0„ + A0«, xn + A Xn +) ), (4.26a) 

with 

A0< +) = -A0<"\ A X « = -A X i") + 2iVvr. (4.266) 

We see from (4.24) and (4.26) that in the rest frame of reference, the asymptotic form 
of the bright iV-soliton solution can be represented by a superposition of N independent 
bright 1-soliton solutions, the only difference being the phase shifts. It follows from (4.24) 
and (4.26) that the formulas for the total phase shifts of the nth soliton are given by 



N 



Ax, r 



E t 

v j=n+l 



Pn + P*j 



Pn - Pj 
24 



n— 1 

E'» 



Pn + P*j 



Pn ~ Pj 



(4.27a) 




(4.276) 

The above formulas reduce to (4.16) and (4.19) for the special case of N = 2. They 
clearly show that each soliton has pairwise interactions with other solitons, namely there 



are no many-particle collisions among solitons. 
5. Concluding remarks 

In this paper, we have presented two different expressions of the bright iV-soliton solution 
of the FL DNLS equation. The exact method of solution developed here is purely algebraic 
which does not recourse to the 1ST. The iV-soliton solution given by theorem 3.1 has been 
obtained by the dressing method [4] while the form presented by theorem 3.2 is new. The 
equivalence of both expressions has been demonstrated by a straightforward computation 
using the properties of the Cauchy matrix. We have also shown that the bright N- 
soliton solution of the DNLS equation can be derived from that of the FL DNLS equation 
through a simple relation. This fact has been stemmed from the analysis of a trilinear 
equation among the tau functions. Last, we have investigated in detail the properties of 
the solutions and found some new features. 

In conclusion, we briefly comment on the dark soliton solutions. The construction 
of soliton solutions under nonvanishing boundary conditions is difficult to perform when 
compared with that under vanishing boundary conditions. Specifically, the analysis by 
means of the 1ST requires the delicate discussion on the spectral problem. On the other 
hand, the bilinear transformation method is applicable easily to obtain soliton solutions. 
In part II of the present study, we will address the construction of the dark iV-soliton 
solution of the FL DNLS equation on a background of a plane wave. The system of 
bilinear equations for the equation is almost the same as that of the bright soliton case. 
However, a difficulty arises due to the constraint imposed on the complex parameters 
Pj (j = 1, 2, N), as in the case of the similar problem for the NLS equation in which pj 
lie on a circle in the complex plane [12]. In particular, this fact must be used explicitly to 
prove the bilinear equation corresponding to equation (2.2). We will show that the trilinear 
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equation analogous to (2.10) play an important role in the proof and the inspection of 
the equation also provides the dark iV-soliton solution of the DNLS equation. 
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